HIERARCHICAL OBJECT FORMATION IN THE PECULIAR 

VELOCITY FIELD^ 

HIDEAKI MOURI 

Meteorological Research Institute, Nagamine, Tsukuha 305-0052, Japan 

o ■ 

^ ■ YOSHIAKI TANIGUCHI 

^ , Astronomical Institute, Graduate School of Science, Tohoku University, Aoba, Sendai 

! 980-8578, Japan; tani@astr.tohoku.ac.jp 

(N 

Cp ■ ABSTRACT 
> 

^ . 

O ' Using the initial peculiar-velocity field, we analytically study the hierarchical 

^ . formation of gravitationally bound objects. The field is smoothed over a scale 

. that corresponds to the mass of a given class of objects. Through the Zel'dovich 

^ , approximation, the smoothed field determines how the objects cluster together to 

'"q_|. form a new class of more massive objects. The standard cosmological parameters 

Q . lead to the evolution of primordial clouds with < 10^ Mq — * galaxies with 10^^ Mq 

0- 



^ clusters of galaxies with 10^^ Mq — * superclusters of galaxies with 10^^ M( 



c§ '. The epochs obtained for the formation of these classes of objects are consistent 

with observations. 

^- 

. Subject headings: cosmology: theory — galaxies: formation — large scale struc- 

ture of universe 



1. INTRODUCTION 

For the formation of gravitationally bound objects such as galaxies and clusters of 
galaxies, the current paradigm is hierarchical formation based on the cold dark matter model. 
The cold dark matter induces significant small-scale fluctuations in the initial field. Since 



^This version corrects errors in equations p9|) . (pO|) . ((22| . and ((25|) and Figures [2}i and [3] of the published 
version (ApJ, 634, 20 [2005]). The errors are also corrected in the erratum (ApJ, 659, 1792 [2007]). 
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the first objects arise from these initial fluctuations, they are not so massive. They merge 
with one another via gravitational clustering in a hierarchy, evolving successively into more 
massive objects. We analyticaUy discuss some statistical properties of this hierarchical object 
formation via the gravitational clustering. 

The existing analytical studies are based on the density fleld, following the pioneering 
works of Press & Schechter (1974; see also Peebles 1980), Bardeen et al. (1986), and Bond et 
al. (1991). They smoothed the initial density field over a scale that corresponds to a given 
class of objects. This is because, even when objects are formed at a scale, the larger-scale 
density field retains its initial pattern, according to the Eulerian linear theory. The smoothed 
initial density field was used to study the gravitational clustering of those objects. 

Our study is based on the peculiar-velocity field, which is expected to offer new infor- 
mation on the gravitational clustering. Over scales where objects have not been formed, the 
peculiar-velocity field retains its initial pattern, according to the Lagrangian linear theory, 
i.e., the Zel'dovich approximation (ZA; Zel'dovich 1970). The initial fiuctuations are random 
Gaussian. For such fluctuations, there are statistical measures (Rice 1954, §3). We apply 
these measures to the smoothed initial peculiar-velocity field. 

The gravitational clustering refiects the cosmological parameters. We use their latest 
standard values: the Hubble constant Hq = 70 km s~^ Mpc~^ {h = 0.7), the cosmological 
constant Q\ = 0.7, the matter density Q^a = 0.3, and the baryon density fl^ = 0.05 (Spergel 
et al. 2003). 

The definitions and equations are summarized in ^ We discuss the general features of 
the hierarchical object formation via the gravitational clustering in §3.1[ The characteristic 
mass scales are discussed in §3.21 The epochs for the formation of objects are discussed in 
§3.31 We conclude with remarks in ^ 

2. DEFINITIONS AND EQUATIONS 

2.1. Eulerian Description and Linear Theory 

The basic assumption for the following equations is that the density field is isotropic 
and homogeneous (e.g., Peebles 1980, §10, 27, and 41). We use comoving coordinates x, 
i.e., physical coordinates divided by the scale factor. The scale factor a{t) is normalized to 
a(to) = 1 at the present time t = to- 

The mass density p{x,t) and its spatial average (p) = 3QraHQ/87iG define the density 
contrast 6{x,t) as p/{p) — 1. Here (■) denotes an average. The relation to the Fourier 
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transform 5{k,t) is 

5{x,t) = j^^^ j 5{k,t)exp{ik-x)dk. (1) 
For the wavenumber k = \k\, the power spectrum P{k,t) is 

{S{k, t)6{k', t)*) = {27rf 6B{k - k')P{k, t). (2) 

Here 6{k,t)* is the complex conjugate of 6{k,t), and ^0(2;) is Dirac's delta function. 

So far as in the linear regime where the linear theory is valid, the density contrast and 
power spectrum grow in a self-similar manner: 

5{x, t) oc D{t) and P{k, t) = D{tfPi^^{k) oc D{tf. (3) 

The linear growth factor D{t) is 



with the Hubble parameter 



g(a) = go(^ + ^^A+ , M . (5) 

Together with equations ([I]) and ([3]), the continuity equation dtp + div(pi;) = yields the 
peculiar velocity 

1 -D(t) ^ik~ 
^(*' ^) = (27r)3/2 y p ^) '^''P^'^ ■ ^^'^^^ 

Here -D(t) is the time derivative of the linear growth factor D(t). The factor kjk'^ implies 
that the peculiar velocity is dominated by larger scales than the density contrast. This 
property of the peculiar velocity is essential to our discussion in ^ 

The two-point correlation of the density contrast (5(33, t)5(a; + r, t)) at the scale r = |r| 
is computed from the power spectrum: 



00 

2 



((5(a;,t)5(a; + r,t)) = 47rD(t)M k' jo{kr) Piin{k)dk. (7) 

Jo 

For the pairwise peculiar velocity, i.e., the relative velocity that is parallel to the separation 
vector r, the two-point correlation {u{x,t)u{x + r,t)) is 



{u{x, t)u{x + r, t)) = ATxb{tf 



3o[kr) - 



kr 



Pnn{k)dk 
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(Monin & Yaglom 1975, §12; Gorski 1988). Here jo{x) and ji{x) are first-kind spherical 
Bessel functions 

• r ^ sin(a;) , . . x sin(x) cos(a;) 

MX) = and ji{x) = . (9) 

We defined the two-point correlations in terms of the linear growth factor D[t) because we 
use them only in the linear regime. 



2.2. Lagrangian Description and ZA 

Consider a mass element that occupies the position at the initial time tin. This 
initial or Lagrangian position a^in is related to the position £c(ajin,t) at the time t as 

Dit) , , 

CC(a;in, t) = CCin + J- -V{Xin, tin), (10) 

according to the ZA (Zel'dovich 1970; see also Shandarin & Zel'dovich 1989). The initial 
velocity 'i;(a3in, tin) is determined on the Eulerian linear theory (eq. [6]). 

Even when the Eulerian linear theory becomes invalid, the ZA is valid (e.g., Yoshisato, 
Matsubara, & Morikawa 1998). Moreover, as in equation (fTOj) . the ZA describes the motion 
of mass elements in gravitational clustering. 

The ZA becomes invalid when the trajectories of mass elements starting from different 
positions intersect one another. The ZA mass elements pass away without any interaction, 
which occurs in practice and leads to formation of a gravitationally bound object (Gurbatov, 
Saichev, & Shandarin 1989; Shandarin & Zel'dovich 1989). We thereby define the object for- 
mation as this trajectory intersection. Since the trajectory intersection proceeds from small 
to large scales (e.g., Yoshisato, Morikawa, & Mouri 2003), the object formation proceeds 
from low to high masses in accordance with the current paradigm (^. The usual definition 
for the object formation is that the density contrast 6{x, t) reaches approximately unity in 
the Eulerian linear theory (e.g.. Press & Schechter 1974; Bardeen et al. 1986; Bond et al. 
1991). This is consistent with our definition because the true density contrast is very large 
in both the cases. 



2.3. Smoothing 



The smoothing is based on a Gaussian window function, which is localized in the real 
and wavenumber spaces and is convenient. If the smoothing scale is R, we equivalently use 
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the new power spectrum 

PUm = exp{-eR'')PUk). (11) 

The relation between the smoothing scale R and the mass Mn of the corresponding class of 
objects depends on the objects themselves (Bardeen et al. 1986; Bond et al. 1991). Since 
we are to assume that the objects have collapsed, we approximate as 

Mn = -R'{p). (12) 

This approximation is in accordance with the definition of characteristic mass scales in §3.21 
(eq. [M]). The mass includes both baryonic and dark matters. 

When we are interested in a class of objects with a mass M^j, the initial peculiar- velocity 
field v{x,tia) is smoothed over the scale R. The smoothed field determines the motion of 
those objects via the ZA (Mann, Heavens, & Peacock 1993). Although the ZA is invalid 
over scales where objects have been formed, it is valid over the larger scales. They are in the 
linear or so-called quasi-linear regime. The ZA on a smoothed field, i.e., the truncated ZA, 
is actually known as a better approximation than the ZA on the original field (e.g.. Coles, 
Melott, & Shandarin 1993; Melott, Pellman, & Shandarin 1994). 

The smoothing over a scale R and the subsequent application of the ZA imply that we 
consider not the initial time but a time when objects with the mass Mr have been formed. 
A larger smoothing scale corresponds to a later time because the object formation proceeds 
from low to high masses (S]). 

The smoothing is also necessary to obtain statistics such as the root-mean-square density 
contrast ( §3.11) . Without the smoothing, in the cold dark matter model, there are significant 
small-scale fluctuations that cause divergence of those statistics (Peebles 1980, §26 and 96; 
Bardeen et al. 1986). 



2.4. Initial Power Spectrum 

The density contrast and peculiar velocity at the initial time are not only isotropic and 
homogeneous but also random Gaussian. Their statistics are uniquely determined by the 
initial power spectrum 

P{k, tin) = D{U,,fPiUk) with PlUk) = T{kfCk. (13) 

The transfer function T{k) is from the cold dark matter model of Bardeen et al. (1986): 

nk) = ^^^^^"^gf^^^ [1 + 3.89g + (16.1g)^ + (5.46g)^ + (6.71g)^] "^^^ . (14) 
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Here the normalized wavenumber q is defined with the shape parameter V in Sugiyama 
(1995): 



9 



Vh 



Mpc ^ with r = r2m/icxp 



2h 



(15) 



We use the scale-invariant Harrison-Zel'dovich spectrum Ck. The normahzation constant C 

is set to be 1 X 10^ on the basis of the root-mean-square mass contrast within the radius 
8h~^ Mpc observed for the present time, 0.9 (Spergel et al. 2003; see also Peebles 1980, §26). 



2.5. Statistics of Random Gaussian Field 

For a homogeneous random Gaussian signal y{x) in a one-dimensional space, there are 
useful statistical measures (Rice 1954, §3). The mean interval between successive zeros is 



A 



71 



{{dy/dxf) 



1/2 



TT 



e(o) 



.1/2 



(16) 



Here ^(r) is the two-point correlation {ii[x)y{x + r)) . The scale A is equivalent to the Taylor 
microscale that is often used to characterize turbulence (e.g., Monin & Yaglom 1975, §12 
and 15). Likewise, the mean interval between successive extrema is 



7] = TT 



\ {{dy/dxf) 1 


1/2 


" d2^(0)/dr2" 


1/2 


_{{d'^y/dx'^Y)_ 


= TT 


d4^(0)/dH_ 
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The probability that y{x) passes through a value y^ with a positive slope in the range 
(x, X + 5x) is 



2A 



exp 



Sx 



2A 



exp 



2^(0), 



Sx. 



(18) 



Since A is the mean size of regions where y{x) has the same sign and there are plus and minus 
signs, the factor exp(— y^/2(y^)) is the probability that \y{x)\ passes thorough a value 
with a positive slope in such a region. This is just the probability that the maximum of \y{x) \ 
exceeds the value We differentiate exp(— y^/2(y^)) to obtain {\y*\/ {y^)) exp{—yl/2{y'^)) 
as the probability that the maximum takes the value 



3. DISCUSSION 



Our discussion is limited to one-dimensional cuts of the three-dimensional peculiar- 
velocity and density fields. Only for such a one-dimensional signal, useful statistical measures 
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exist ( §2.51 and 13.21) . They have been nevertheless used in three-dimensional turbulence for 
decades (e.g., Monin & Yaglom 1975, §12 and 15) and are hence expected to yield crude but 
still useful information on three-dimensional process of the object formation. 

3.1. General Features 

Figure [1] shows the peculiar velocity Vx{R\x,ti^) and the density contrast 6{R\x,ti^) 
at the initial time tin on a one-dimensional cut x smoothed over the scale R = 0.01 Mpc. 
They are normalized with their root-mean-square values. The peculiar velocity Vx is in the 
direction of the one-dimensional cut. 

The peculiar velocity and density contrast were obtained from the three-dimensional 
Fourier transforms. Their real and imaginary parts are mutually independent and random 
Gaussian. We assigned the Fourier transforms according to the initial power spectrum (eq. 
[2]) and converted them to the real-space counterparts (eqs. [T] and [B]). 

The density field has many peaks. They individually collapse into gravitationally bound 
objects. The peculiar- velocity field has larger-scale variations ( §2.ip . Then, as predicted by 
the ZA and indicated by arrows, the density peaks cluster together. They form a new class 
of more massive objects, provided that the clustering velocity is high enough]^ Zeros of the 
peculiar velocity, i.e., stationary points of the clustering motion, are centers or edges of the 
clustering. There the density contrast is positively or negatively enhanced in a significant 
manner. We also see small-scale velocity variations, which correspond to collapse and merger 
of the individual density peaks during the clustering. 

While the mean size of density peaks is equal to the mean interval rjs between successive 
extrema of the density contrast, the mean size of overdense regions is equal to the mean 
interval Xs between successive zeros of the density contrast. The mean clustering scale is 
equal to the mean interval between successive ridges and troughs of the peculiar velocity, 
where the clustering motion is largest 1^1 This interval is equal to the mean interval A„ between 
successive zeros of the peculiar velocity. Together with the typical clustering velocity, the 
mean clustering scale determines the typical clustering timescale. 



^ The collapse of the individual density peaks into gravitationally bound objects occurs earlier than the 
clustering into a new class of more massive objects. 

This interval does not represent edge-to-edge sizes of the clustering regions. We thus ignore their outer 
envelopes. 
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3.2. Characteristic Mass Scales 



Using equation ([8]) and its derivatives at r = 0, the length scales Xu{R) and rju{R) in 
equations f|T6|) and f|T71) are obtained for the pairwise peculiar velocity at the initial time tin 
smoothed over the scale R: 



and 



Xu{R) 



Vu{R) 



71 



IT 



5 j;;^ PiUR\k)dk 

k^PrUR\k)dk 

7 k^PiUR\k)dk- 
.5 k'^Pii^{R\k)dk^ 



1/2 



1/2 



(19) 



(20) 



Likewise, using equation ([7]) and its derivatives at r = 0, the length scales Xs{R) and r]s{R) 
are obtained for the initial density contrast: 



and 



XsiR) 



vs{R) 



3 ePun{R\k)dk 

j^~k^K(mdk^ 

5 k^Pnn{R\k)dk' 
.3 r k^Pu^{R\k)dk^ 



1/2 



1/2 



(21) 



(22) 



The scales riu{R) and XsiR) are identical apart from a numerical factor. For the initial density 
contrast, we also obtain the correlation length ks{R), i.e., the length scale characteristic of 
the large-scale correlation (e.g., Monin & Yaglom 1975, §12) 



tis{R) 



j;^{6{R\x,U^)6{R\x + r,U^))dr _ vr kP^^{R\k)dk 



{5{R\x,U^Y) 



2r k^PrUR\k)dk' 



(23) 



The correlation length is zero for the initial pairwise peculiar velocity because its two-point 
correlation is significantly negative at large scales (see below). Although Xu{R), Vu{R), 
Xs{R), ris{R), and ns{R) are defined at the initial time, they are valid throughout the linear 
regime where the peculiar velocity and density contrast remain random Gaussian. 

The mass scales Ma, M^, and are defined as the masses of spheres with the diameters 
A, T], and n: 

3 



An 



(p)- 



(24) 



The mass scales Mx^(^r) and Mr^^i^R) are Lagrangian quantities and are hence valid so far as 
the ZA is valid. We regard the mean clustering mass Mx^(r) as the mass for a new class of 
objects that are formed by objects with the mass Mr. 
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Figure [2]2 shows the mass scales for the pairwise pecuhar velocity Mx^(r) and M^^(^ji) 
and those for the density contrast Ma^(r), M^^(jj), and M^^^f^n) as a function of the smoothing 
mass Mr. The corresponding length scales are also shown. 

Since the mean size of density peaks ris{R) is comparable to the smoothing scale R, the 
density peaks are the objects that correspond to the smoothing (Bardeen et al. 1986)0 The 
mean size of overdense regions Xs{R) is a few times greater. Since the mean clustering scale 
Xu{R) is far greater than the mean size of density peaks rjslR), the density peaks cluster 
together to form a new class of more massive objects. Since the mean clustering scale Xu{R) 
is also far greater than the density correlation length ks{R), the clustering is independent of 
the individual density peaks. 

The clustering is not significant for the length scales A„(i?) > 10^ Mpc, where we see 
Xu{R) — Xs{R). While the individual overdense regions do not cluster together but collapse 
to single objects, the individual underdense regions expand to single voids. They form the so- 
called cellular structure where the objects constitute the cell walls. The two-point correlation 
of the pairwise pecuhar velocity {u{x,ti^)u{x + r,ti,^)) at r > 10^ Mpc is correspondingly 
negative as in Figure 

Figure [2I2 is used to discuss the hierarchical formation of gravitationally bound objects, 
via successive smoothing of the peculiar-velocity field. The first objects are clouds with 
10® Mq, which are considered to harbor the first stars in the universe (e.g., Bromm & Larson 
2004). The smoothing for Mr = 10® Mq yields the mean clustering mass M\ ~ 10^^ M©, 
i.e., the typical mass for galaxies (e.g., Brainerd, Blandford, & Small 1996) 1 Thus primor- 
dial clouds form galaxies. The smoothing for Mr = 10^^ Mq yields the mean clustering 
mass Mx^(R) — 10^^ Mq, i.e., the typical mass for clusters of galaxies (e.g., Merritt 1987). 
Thus galaxies form clusters of galaxies. The smoothing for Mr = 10^^ Mq yields the mean 
clustering mass Ma„(r) ~ 10^® M©. Thus clusters of galaxies form superclusters of galaxies. 
They are practically the most massive gravitationally bound objects because the cluster- 
ing timescale for Mx^^r) > 10^® Mq is too large f rojl . The mass scale Mx^(r) ~ 10^® M© 
corresponds to the length scale Xu{R) — 10^ Mpc. 

The smoothing mass Mr in our above discussion was set to be the typical mass for an 
actual class of objects. Thus our discussion was analogous to observations, where stars are 



The scale rju{R), which is considered to represent the motion associated with the coUapse of the indi- 
vidual density peaks, is also comparable to the smoothing scale R. 

^ The mean clustering mass Mx^[r) is insensitive to the smoothing mass Mji if it is less than about 
10^ Mq. This is because, as seen in equation p^ . the mean clustering scale A„(i?) is not so dependent on 
small-scale fluctuations. 
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used to identify galaxies and galaxies are used to identify clusters of galaxies. It should be 
also noted that we considered a time when the class of objects with the mass have been 
formed (§23]) • 

Our success in reproducing distinct classes of gravitationally bound objects, i.e., galax- 
ies, clusters of galaxies, and superclusters of galaxies, was due to the existence of the mean 
clustering scale \u{R) or the mean clustering mass Mx^(r) that is distinctly greater than the 
smoothing scale R or the smoothing mass Mr, respectively. We underline that Xu{R) and 
Mx^(R) are from the peculiar- velocity field. The density field does not offer information as to 
how the objects cluster together to form a new distinct class of more massive objects, at least 
at the level of the standard Eulerian linear theory. This is because, although the density field 
is related to the peculiar-velocity field, it is dominated by smaller-scale fluctuations (Fig. [H 
see also §2.11) . 



3.3. Formation Rate of New Objects 

Within our framework of hierarchical object formation ( §3.2p . we estimate the time for 
objects with a mass Mr to form a new class of more massive objects. The mean clustering 
scale is A„(i?). The mean clustering mass, i.e., the mass for the new class of objects, is 

The initial clustering velocity of objects with a mass Mr is of the order of {u{R\x^ tin)^)^"^^. 
According to the ZA (eq. [ID]), the displacement of the objects by the time t is D{t)/D{ti^) 
times 

i/(M(i?|x,tin)2)V2. Here z/ is a nondimensional parameter. When the trajectories of the 
objects intersect one another, the objects cluster together to form new objects ( §2.21) . Since 
the mean clustering scale is \u{R), the typical radius of the clustering region is \u{R)/2, 
which is regarded as the typical displacement. Using equations ([8]), f|T6l) . and f|T9l) . we have 



m = ^ 



57r 



J^'^k^PyUR\k)dk\ 



1/2 

(25) 



For a value of the parameter p, equation (l25l) determines the formation time of the new 
objects with the mass Mx^(r). 

Figure [3h shows the mass Mx^{r) of new objects formed at redshift z(t) + 1 = a(t)~^ 
for the parameters u = 1-5. The mass is high at a low redshift. At each redshift, the mass 
is high for a large u value, i.e., a large clustering velocity, which tends rare. Objects with 
Mx^(R) ^ 10^^ Mq are not formed because the object formation is suppressed at z < 0. 
This suppression is due to the accelerated expansion of the universe. The increase of the 
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linear growth factor D(t) is no longer significant at z < as indicated by the dotted curve 
in Figure©. 

Given the mass M^, equation fl2^ determines the parameter as a function of the time 
t for the formation of new objects with the mass Mx^{r)- The parameter u represents the 
maximum velocity in the clustering region. Then its probability p{v)5v is i/exp(— z/^/2)5z/ 
( §2.5p . Using the relation 5v/v = 6D/D, which is from equation ( !25l) . the probability is 
rewritten as 



Equations (12^ and fl2B]) determine the formation rate of the new objects. 

Figure [3l> shows the formation rates as a function of redshift z{t) + 1. For galaxies with 
MA„(ij) ~ 10^^ Mq, their formation by primordial clouds with M/j = 10^ Mq starts at 2; ~ 10 
and is maximal at z ~ 3. For clusters of galaxies with Mx^^r) — 10^^ Mq, their formation 
by galaxies with Mr = 10^^ Mq starts at ^ ~ 2 and is maximal at 2; ~ 0. This maximum 
corresponds to the saturation of the linear growth factor D{t). For superclusters of galaxies 
with Mx^{R) — 10^^ Mq, their formation by clusters with Mr = 10^^ Mq starts at z < 0. 

The formation rates in Figure [3^ are consistent with observations. Since the star for- 
mation rate integrated over galaxies has a broad peak at 2; ^ 4-1 and declines to its present 
value (e.g., Heavens et al. 2004; Taniguchi et al. 2005), the typical galaxy mass is assembled 
by 2; ~ 1. The number density of clusters of galaxies with > 10^^ increases by a factor 
of 10 from 2; ~ 1 to 2; ~ (Bahcall & Bode 2003). Since the present galaxy distribution is 
almost uniform over scales above 40Mpc (e.g., Kurokawa, Morikawa, & Mouri 2001; Jones 
et al. 2004), superclusters of galaxies have not undergone significant gravitational clustering 
at z ~ 0. 

There are analytical studies of the formation rate based on the density field (e.g., Lacey 
& Cole 1993; Kitayama & Suto 1996). The formation rate was obtained by integrating the 
transition rate from a mass M to another M'. Since the transition rate diverges in the limit 
M — > M', the formation of an object was defined as the assembly of a half of its mass. 
This definition is artificial and also does not correspond to the situation that objects cluster 
together to form a new distinct class of more massive objects. The density field does not 
offer the required information ( §3.2p . We prefer our approach based on the peculiar- velocity 




(26) 



field. 
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4. CONCLUDING REMARKS 

The initial peculiar-velocity field was smoothed over a scale R that corresponds to the 
mass Mr of a class of gravitationally bound objects. Through the ZA, the smoothed field 
determines how the objects cluster together to form a new class of more massive objects 
(Fig. [T]). The mean clustering scale Xu{R) yields the mean clustering mass Mx^(r), which is 
regarded as the mass for the new class of objects. We assumed the standard cosmological 
parameters and obtained the evolution of primordial clouds with < 10^ M© — > galaxies with 
10"*^^ Mq — > clusters of galaxies with 10"*^^ Mq — > superclusters of galaxies with 10^^ Mq (Fig. 
I2b). Together with the clustering velocity, the mean clustering scale A„(i?) was also used to 
obtain the epochs for the formation of those classes of objects (Fig. [2t). They are consistent 
with observations. 

The peculiar- velocity field allowed us for the first time to analytically study the forma- 
tion of distinct classes of gravitationally bound objects. This is because the peculiar-velocity 
field is dominated by large-scale fluctuations (Fig. [T]). There is the mean clustering scale 
Xu{R) that is distinctly greater than the smoothing scale R. Since it is difficult to obtain such 
a characteristic from the density field, the existing analytical studies based on the density 
field did not reproduce the distinct classes. There is even a model of continuous clustering, 
where galaxies and superclusters of galaxies are extremes of a continuous range of objects 
(Peebles 1980, §4). This model does not appear so reahstic. We also underhne that the ZA 
used in our study is a better approximation than the Eulerian linear theory used in most of 
the existing studies. 

Our study is based on one-dimensional cuts of the three-dimensional peculiar-velocity 
field. Although Bardeen et al. (1986) obtained useful statistical measures for a three- 
dimensional density field, there is no useful statistical measures for a three-dimensional 
velocity field. The advent of such measures is desirable. Then we could obtain more infor- 
mation on the peculiar-velocity field and hence the hierarchical formation of gravitationally 
bound objects. 

The authors are grateful to M. Morikawa and A. Yoshisato for interesting discussion 
and to the referee for helpful comments. 
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Fig. 1. — Peculiar velocity Vx{R\x,tin) and density contrast S{R\x,tin} at the initial time tin 
on a one- dimensional cut x smoothed over the scale R = 0.01 Mpc. They are normalized with 
their root-mean-square values. The velocity is in the direction of the one- dimensional cut. 
The arrows indicate the clustering motion. We also indicate an example of the clustering 
scale. 
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Fig. 2. — (a) Mass scales for the pairwise peculiar velocity Mx^{r) and M^^(ij) as well as mass 
scales for the density contrast MA^(i?), Mrig(R), and M^^i^r) as a function of the smoothing 
mass Mr. The corresponding length scales arc also shown. (6) Two-point correlation of the 
pairwise peculiar velocity (^(a;, iin)M(£C + r, tin)) normalized with (^(a;, ^m)^). 
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Fig. 3. — (a) Mass Mx„{r) of objects formed at redshift z{t) + 1 for the clustering-velocity 
parameter u — 1-5. The length scale Xu{R) and the scale factor a{t) are also shown, (b) 
Formation rates at redshift z{t) + 1. We consider galaxies with Mx^(r) — 10^^ Mq formed by 
primordial clouds with Mr = 10^ Mq, clusters of galaxies with MA„(i?) ~ 10^^ Mq formed by 
galaxies with Mji = 10^^ Mq, and superclusters of galaxies with Mx^{r) — 10^^ Mq formed 
by clusters with Mr = 10^^ Mq. The dotted curve denotes the linear growth factor D{t). 



